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^ \ We consider problems of one dimensional interacting fermions confined to a finite size, multichan- 

Q . nel geometry. Concentrating on Luttinger liquids and carbon nanotubes, we use nontrivial boundary 

' conditions to represent the effect of external leads, and apply our framework to transport problems 

' in a Josephson junction setup. We present an exact computation of all correlation functions, includ- 

fT^ , ing finite-size and temperature effects, for two sets of solvable boundary conditions. In all cases, we 

compute physical quantities like the Josephson current and the pairing order parameter profile. 

(D ■ I. INTRODUCTION 

u . 

c/3 ' It is now common knowledge that one-dimensional interacting electronic systems possess metallic phases that 
are fundamentally different from those of their higher-dimensional brethren, which are for most practical purposes 

■ extremely well described by adaptations and refinements of Fermi liquid theory. The Luttinger liquid has become the 
paradigm for the former systems: the quasiparticle pole is no more, and instead one finds only collective spin and 
charge excitations living around two Fermi points. The most striking consequence of interactions are first that spin 
and charge excitations are not bound to one another anymore, and can thus travel at different velocities. Moreover, 

Q . charge fractionalization can occur, whereby "fundamental" quantities like the electron's charge can be split in many 
Q \ pieces (for a review, see e.g. 

Experimental realizations of Luttinger liquids range from edge states in the fractional quantum Hall effect (in 

^ ^ this case, a chiral Luttinger liquid is found), quantum wires in semiconductor heterostructures and perhaps most 

^ , notably in sing le- walled carbon nanotubes (SWNT) ||, |, |. The latter can be used to perform electrical transport 

■ experiments using different geometries, for example in junctions manipulated by mechanical means |@, |^, ^ . "Kinks" in 
' the nanotubes or crossings of nanotubes can be produced by scanning tunneling microscopes (STM) tips, and these in 
> turn provide realizations of backscattering impurities embedded in a Luttinger liquid or of connections through tunnel 
' junctions. These junctions can act as rectifying diodes at room temperature, thus opening up a route towards the 
. construction of nanoscale devices of all sorts. An example of recently proposed applications is that of an "entangler" , 
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whereby two nanotubes coupled to a superconductor are expected to produce physically separated entangled pairs of 
excitations 0, |ll|, |l|| . 

Single impurities embedded in Luttinger liquids have attracted a lot of theoretical interest since the early work of 
Kane and Fisher [ p^ . In particular, many app roaches have been devised based on boundary conformal field theory 
I (CFT) iQ and boundary integrability Q. The trick here is to "fold" the system around the impurity, thereby 
■ O replacing it with an effective boundary whose effects scale either to the weak- or strong-coupling regime depending 
^ ; on the nature of the bulk interactions in the system. For example, for a backscattering impurity in a Luttinger liquid, 
repulsive interactions enhance the backscattering amplitude, and therefore suppress the conductance. 

Another extremely interesting situation to consider is that of a Luttinger liquid on a finite interval. This class of 
problems is a form of generalization of the previous single impurity ones, in the sense that we now have two different 
boundaries to deal with. One thus imagines electrons being trapped by either a direct physical cut of their support 
(maybe an ended nanotube), or prevented from propagating further by the presence of an excitation gap, maybe a 
\ superconducting gap, or voltage gates. We consider here the problem of a Josephson junction realized by bridging two 
superconductors with either quantum wires or carbon nanotubes. Experimental attempts at realizing this device can 
be found in Theoretically, this problem has been addressed in 18 1, where the Josephson current for "perfect" 
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contacts was computed (here, by perfect, we mean a contact such that only Andreev reflection occurs; a poor contact 
would be one such that most of the reflection processes that occur would be normal reflection ones). A perturbative 
analysis of the problem was presented in [|l9|, and an effective model amenable to boundary CFT was used in po|. 
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For certain cases, the problem is an integrable quantum field theory, and this was treated in However, only 

one particular (nontrivial) value of the interaction parameter can be fully solved, the others leading to rather tedious 
complications with the analytical properties of the thermodynamic Bethe Ansatz, which are yet to be fully resolved. 
One common factor of all these studies, however, is that they were devoted to single-channel Luttinger liquids. 

In the case of nanotubcs, the number of channels available for conduction is however by construction greater than 
one. Moreover, if one wants to describe interesting devices like the entangler mentioned above, or maybe even things 
like multiwall carbon nanotubes [ p2[ , it is important to have a theory for multichannel situations. It is very tempting 
to think that the one-dimensionality of the system opens the way to a full nonperturbative solution of transport 
problems of this class, as it can do for the single-channel case. However, from basic factorized scattering arguments, 
it is possible to show [ p3| that multichannel theories nontrivially coupled only at boundaries but not in the bulk, do 
not seem to have enough flexibility to be integrable (that is, only a trivial scattering matrix can satisfy the bootstrap 
equations). This holds, of course, if one considers having no boundary degrees of freedom. Relaxing this might lead 
to some other integrable theories, but not ones which can be directly applied to the problem at hand. One is therefore 
ultimately led straight back into the domain of perturbative computations as the only fully controllable approach, 
and most of what we will present here shall be based on a minimal use of it (that is, we solve all that we can exactly, 
using perturbation theory only when strictly needed). 

Our plan is therefore to treat the general problem of a conducting multichannel system in a finite size domain, 
between two boundaries having nontrivial effects on the physics (but no independent degrees of freedom of their own). 
The noninteracting and interacting cases are treated in fundamentally different manners. The noninteracting problem 
is straightforward to solve thoroughly directly in the fcrmionic language, providing in some sense an adaptation of 
the well-known work of Blonder, Tinkham and Klapwijk to two-boundary problems. The partition function can 
be computed exactly, which means that thermodynamic properties and correlation functions can be obtained exactly 
(for the case of correlation functions, we only solve the interacting case, since the noninteracting ones can be recovered 
easily by putting all interaction parameters to zero). 

Including interactions, however, takes us to a different set of problems. First of all, the presence of appropriate 
Coulomb interactions renders the fcrmionic description untractable. Fortunately, however, we can use bosonization 
to perform our computations. This, as is well-known, allows one to treat interactions nonperturbatively. In the bulk, 
correlators are rather easily computed. Remember, however, that we are here dealing with two-boundary problems. 
This involves treating Luttinger liquids and carbon nanotubes in a finite size at finite temperatures, for which we 
present a computation of all multipoint correlation functions around simple fixed points represented by tractable 
boundary conditions. The scaling of the correlators then depends on whether the operators sit away from or near 
to the boundaries: the scaling near the boundaries differs from that in the bulk, and depends on the specifics of the 
boundary conditions that are in use. We venture to treat this problem in all generality at the level of bosonization. 
For our specific problems, we set up the perturbative formalism, and present computations for physical transport 
quantities like the Josephson current in different situations, as well as discussions about the pairing order parameter 
profile. 

The plan of the paper is as follows. In section 2, we treat the noninteracting problem in fcrmionic language, and 
compute finite-temperature transport properties using canonical mode expansion methods. In section 3, we consider 
the interacting case. We consider separately the cases of boundary-coupled Luttinger liquids, and nanotubes. For each 
of these, we discuss separately the cases of normal and Andreev boundary conditions, representing "bad" and "good" 
coupling to the superconductors (again, in the sense described above in terms of normal versus Andreev reflection 
amplitudes; the microscopies of this is a rather tricky experimental issue). The physical results are listed in this part 
of the paper, and summarized in the conclusion. Important formulas for the bosonization, including computations of 
correlators at finite temperature in a finite-size geometry are collected in Appendix A. In Appendix B, we provide 
some definitions of (multivariable) 0-functions appearing in our formulas for the correlation functions. Appendix C is 
devoted to some notes on the real-time, finite-temperature perturbative formalism we used. 

II. NONINTERACTING CASE 

In this section, we concentrate on the case where interactions are absent in the bulk of the system. The main 
advantage of this is that the calculation of thermodynamic quantities and of correlators can be performed exactly 
using relatively straightforward methods. We choose to always work in a real-time formalism, even in the presence 
of finite temperatures. This allows an immediate adaptation of our framework to time-dependent perturbations 
(although we do not consider them here, they are in principle easy to address using our formulas, bearing the notes 
in Appendix C in mind). Moreover, the definitions and results obtained in this section, in particular the treatment of 
boundary phenomena, will prove to be a good starting point for other problems later on, when we will consider other 
physical cases with interactions in the bulk. Our intention is to provide a flexible framework for building theories of 
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e.g. transport through one-dimensional channels, in many different types of setups. Recent advances in experimental 
realizations of nanostructures provide good motivation for such types of theories. 

Our starting point is a theory for A'' decoupled spinful channels of fermions living within a one-dimensional channel 
of finite size, x £ [0, i?], which we often call the fundamental domain. The spin-1/2 Fermi fields are labeled by channel 
i = 1,...,N and spin indices, 'I'^(a;). The canonical equal-time anticommutation relations for the fermions are (for 
x,x' e]0,R[) 

{vE'i(x,t),<;(x',t)} = S^^,P^'S{x-x'). (1) 

In order to write proper mode expansions for our fermions, we have to be careful with what happens in the presence 
of the two boundaries. The crucial thing to preserve is the independence of ^'^(a;) at both boundaries: we do not 
want to identify fields at a: = and fields at x = R, as one would normally do using a finite-size setup on a ring. 
We here have to start by unfolding the fields, and then requiring some form of periodicity. The way we do this is by 
extending the definitions of the fields to the whole real axis by reflecting the original field at both ends of the system, 
i.e. by taking 

^i{~x)^^i{x), ^iiR-x)^^i{R + x). (2) 

This immediately shows that the proper periodicity of the fields is 2i?, and not R as we would have obtained in a 
closed ring of circumference R. 

Each channel has a value of the Fermi wavevector kp and Fermi velocity which we take as external parameters. 
This gives us the flexibility for example to treat the cases where the different channels are put at different external 
voltages In the simplest case, the relationship between the external voltage applied to channel j and the 

Fermi parameters is kp = kp + where kp is the Fermi wavevector at half- filling. This step in not crucial at this 
stage, as most of our formulas will comprise boundary backscattering potentials representing the voltages of leads, so 
we can simply keep the voltages explicitly instead of reabsorbing them. For transport, only voltage differences matter. 
The Fermi fields are then linearized as usual in terms of chiral left- and right-movers as 

- e^'-^^'j,,{x) + e-^''>^i^{x). (3) 
with canonical equal-time anticommutation relations 

{^Ux,t),^il,{x\t)} = (4) 

and a similar equation for the right movers (the left and right fields anticommute in the fundamental domain; note 
the factor of a half, which is usual for chiral fields |^). In terms of left- and right-movers, the extensions (^) imply 
the analytical continuations of chiral fields 

n.(.-^) = KA^): nAR-^) = e'^'^'^^^^^jR + x), (5) 

which represent a full mapping between the left- and right-movers in the theory (that is, we can use one or the other, 
as long as we choose the right periodicity for getting the proper mode expansions). The appropriate periodicity of 
the fields then becomes 

^i^{2R + x)^ e^.fe^-^,*^^^^^^ KA^R + ^)= e-'''-^*fl.(x). (6) 

The above formulas in turn imply what we term the normal boundary conditions at the left and right ends of the 
system, representing the fact that excitations hitting the boundary are simply normal-reflected back into the bulk of 
the system: 

*i.(0) = *L(0), = e^^'^'^'^^^^JR). (7) 

The important point to realize is that these boundary conditions represent a link between the fields of the theory: 
using them, all the fields of the theory can be described either by left- and right-movers in the fundamental domain, 
or using only e.g. left-movers, defined over an interval twice as big as the fundamental domain. This set of boundary 
conditions is used throughout (half) of the paper, and represents a basis from which bosonization can be done exactly. 
They should be kept closely in mind by the reader. 
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The action describing the left- and right-movers with normal boundary conditions is composed of a noninteracting 
bulk part and a specific boundary part which is chosen such as to canonically impose the normal boundary conditions 

(0): 



CI q{bulk) q(bdry) 

^norm — '->Q ~r '->noT"m 



(8) 



The noninteracting bulk action is most conveniently written in terms of a symmetric form of chiral left- and right- 
movers: 



S 



(bulk) 



R N 

dt I ^ 







(9) 



The boundary action imposing normal boundary conditions is 



(10) 



This theory is very simple: its partition function and correlators are straightforward to obtain, even at finite temper- 
ature. One can easily verify that the equations of motion at the boundaries coincide with the boundary conditions 
(0)- 

What really interests us, however, is to couple the system to external leads. One of the most interesting configu- 
rations that we can consider is the one where the system is in a Josephson junction geometry, whereby the left and 
right ends are in contact with two superconductors having different order parameters. On general grounds, one then 
expects currents to flow through the system even in the absence of voltage bias, part of the phenomena called the 
Josephson effects . 

We therefore consider coupling our system to two external superconductors with bulk gaps well above the excitation 
energies existing within the conducting channel. In particular, this means that quasiparticle penetration within the 
superconductors from the channel are basically non-existent, and that the excitations are either normal- or Andreev- 
refiected immediately back into the channel by the contacts with the superconductors. An effective theory is therefore 
obtained by integrating out the superconductors themselves po| , leaving behind effective boundary actions taking 
the form of BCS-like couplings at the ends of the system. Therefore, we perturb the free action by the additional 
boundary actions 



S 



Abdry) 



(bdry) 



N 



CT 



(11) 



The boundary BCS couplings A^ jj are effective couplings coming from the integrating out procedure detailed above 
(we refer the reader to the discussion in paper pO| , where this is presented in all necessary details) . It is very important 
not to confuse them with the bulk couplings of the original superconductors: the boundary pairings contain information 
about the quality of the contact, the density of states of the superconductor at the junction, etc. We take them as 
external parameters. Their phase, however, can be reasonably supposed to be that of the original superconductors at 
low energies. 

The one advantage of this form of the boundary perturbation is that it preserves the bilinear form of the action as 
far as the fermions are concerned. The theory can therefore again be solved exactly. Varying the action yields the 
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interpolating boundary conditions 



J o" 



-ik] 



(12) 



2 o"' 



where e^^ = — ej,| = 1- These boundary conditions aUow us to move continuously between various extreme cases, 
where things simplify considerably. Besides the normal boundary conditions above, we also can tunc our boundary 
parameters to achieve what we will call Andreev boundary conditions: taking 



= 0, 

the Andreev boundary conditions read explicitly 

-,Jt 



(13) 



(14) 



The two cases of normal and Andreev boundary conditions are those for which one can compute thermodynamic 
quantities exactly in the interacting case, and we urge the reader to keep them in mind in what follows. Away from 
these, one has to rely on some form of perturbation theory. 

In the absence of interactions, however, the exact solution of the theory is obtained by simple considerations. Since 
the fermions are free chiral fields in the bulk, we can write the mode expansions 



Substituting ( p^ ) in (|2|) and defining the Nambu spinors 



(15) 



V 



(16) 



we get the continuity equations 



[1 - iv-^Mi] = [1 + iv-^Mi] ^^fRiuj), 



where v — diag(ti'), and the boundary matrices are given by 



,i<Tj(gikF-R+ili2iv '-ujR 



M 



L,R 



yL,R Al,R 

A* _V* 
^L,R ^ L.R 



(17) 



(18) 



(note: we use here a Bogoliubov-de Gennes notation, with each entry of the 2x2 matrices above being an NxN 
matrix, where N is the number of channels. We order the tensor products as (BdG space) ® (channels space)). The 
quantization condition is then obtained by the requirement that a given eigenstate energy uj fulfills both continuity 
equations. This implies 



det [1 + F{uj)] = 0, 



with matrix 



F{uj) ^ -[l-i-v-^ML] ^[l + iv-^Mije 
X [1 - Iw-^Mr] [1 + iv^^Mfl] 



(19) 



(20) 
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Noting that we can write 

[l-zv-^M]"^ [l + zv-^M] = e^'^^^^Wv-^M)^ (21) 

we see that = F~^, so all solutions to the quantization equation are such that w G M. 

The ground-state energy of the system is then given simply by the sum over a Fermi sea, which we can write as a 
contour integral: 

Eo = y^u) = / divco , , = / dwlndet[l + F(lu)] (22) 

<^o ^""'Jc det[l+F(a;)] 27Ti Jc ^ V ;j V ; 

where the contour C circles all the poles of det[l + F{oj)]~^ on the negative real axis, in a clockwise direction. Using 
complex plane manipulations, we can write 



Ec) = — 
" 27ri 



-1 ptOO -I pOO 

— dLulndet[l + F{uj)]=-— dr/lndet ([1 + i^(-zr?)][l + F(ir;)]) = 
7r« J-ioo 27r Jo 

1 r°° 

^~2^Jo ^^liidet([l + F(-ir?)][l + F-i(ir?)]) + (23) 

where the first term is now manifestly convergent, and the additional contributions are the boundary intensive energies, 
which we ignore from now on since they have no influence on the physics that we are after. 

We can now introduce finite temperatures by inserting temperature-dependent occupation numbers for all states, 
and recalculating the sum. We find, for the equilibrium distribution, 



InZ 



/OO 
dojD{oj)ln 1 + - l3Eo (24) 

-oo 



where the thermal density of states is given by the expression 

1 d 

D((j) = - lim — In det [1 + F((j- iS)] (25) 

where denotes the imaginary part. A more useful expression for numerical purposes is obtained by closing the 
contour around the poles of the Fermi distribution. We then get 

oo 

lnZ = ^Trln[[l + F(-27r(n + l/2)i//3)][l + ^^-i(27r(n + l/2)V/3)]] • (26) 

Although the sum cannot be taken in closed form for general boundary couplings, higher-order terms become expo- 
nentially suppressed at a given finite temperature. For zero temperature, the sum becomes an integral reproducing 
the formula above. 

The crucial point to bear in mind here is that we now have an expression for the partition function, for finite size 
and finite temperatures, for an arbitrary set of boundary parameters. These, hidden in the matrices above (which 
contain the boundary matrices), allow one in principle to compute many thermodynamic quantities in the system by 
using simple parameter derivatives. We will in what follows mostly concentrate on a particular one, but the reader in 
urged to remember that the formulas above are more adaptable than we shall make out. 

Let us now turn to the computation of interesting observables, the most important of which is the total current 
through the system. This is given by the expectation value of the operator (we put the electron's charge to one) 



It{x) = iJ2y' [*L(a^)*i^o(^) - *L(a^)*«-H ■ (27) 



jo- 
in terms of chiral components. Now suppose that for all ij, we have A^l = |A^|,A^ = |A^|e~''^. The physical 
interpretation of this is clear: the coherence length of the original superconductor is large enough to subject all 
channels to the same pairing phase {i.e., we have seen that the bulk pairing phase gets transferred more or less 
directly to the phase of the boundary pairing. If all the channels are coupled to the superconductor within a domain 
of typical size smaller than the coherence length, then the induced boundary pairing phase is the same for all channels). 
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12 3 

FIG. 1: The Josephson current at zero temperature for different values of the boundary parameters as a function of the phase x- 
This is for the 2-channel case, with no interactions, and with equal pairing amplitudes (that is, all the pairing matrix elements 
are equal). The critical current at zero temperature in this figure and the one below has been put to 1. 



The modulus of the boundary pairings, however, can still vary from one channel to the other: microscopic details 
then have influence over them. In that case, it is straightforward to show the well-known fact that 



ill 



2 d 

7dx 



■InZ 



(28) 



using sources and an x-dependent gauge transformation. Therefore, a knowledge of the phase dependence of the 
partition function allows one to compute the current without further work, simply by using this derivative trick. 
This derivation is valid even in the presence of interactions, as long as these can be expressed as current-current 
interactions. 

At zero temperature, we obtain the simple formula 



Ux) 



d 

2—Eo 

ax 



(29) 



involving the ground-state energy above. This is plotted in figure (^ for the case of boundary pairing with equal 
entries, i.e. A'^ — A. This essentially reproduces the plot in EG], done for a single channel, the only difference being 
that the perfect current occurs for A = 1/2 instead of one. One can clearly see how the shape of the current-phase 
relationship changes from a sinx behaviour for small boundary pairing (so with normal reflection amplitude >> 
Andreev reflection amplitude) to the sawtooth function x (mod 2tt) for optimal boundary pairing (where Andreev 
reflection amplitude >> normal reflection amplitude). 

At finite temperatures, we can obtain a closed form expression by computing the ^-derivative by hand. Defining 
the functions 



we can then write 



arctan(v 



'Ml) -icr^igikFfl-lgiv^irji? • 



^2iarctan(v M_n.) 



oo 

Hx) = [[1 + F{-if]n)]-'F{-i7^n) + h.c. 



(30) 



(31) 
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FIG. 2: The Josephson current at different temperatures for a few values of the boundary parameters as a function of the phase 



where r]n = ^{n+ 1/2). Such a current is plotted in figure (|^) for different temperatures and for a few values of the 
boundary parameters. 

The appearance of the commutator in equation (^) is easily interpreted: the terms in (|2^) that do not commute 
with the (Tz matrix are precisely the boundary pairing terms driving Andreev reflection, and therefore the Josephson 
current. 

This completes our discussion of the noninteracting case. From our formulas, it is easy to compute other transport 
properties, or to adapt to different systems with, for example, boundary terms breaking the spin symmetry. We leave 
these for now, and turn instead to the effects of interactions in the bulk. 



III. THE INTERACTING CASE 



In this section, we study two general classes of interacting one-dimensional systems at finite size, coupled to the 
boundary conditions that simulate the superconducting leads to which the real system is coupled. We consider first 
the case of a multichannel Luttinger liquid in which the different channels only interact through the boundaries. Then 
we study the case of a carbon nanotube which can be modelled as a two channel Luttinger liquid with a particular 
choice of interchannel coupling. In both cases, we consider as the starting point the system coupled to either normal or 
Andreev boundary conditions, and introduce more general boundary conditions perturbatively. These lead to different 
types of behaviour, which we detail by presenting computations of the Josephson current (either perturbatively or 
exactly depending on the situation), and of the superconducting order parameter profile in the bulk of the system 
induced by the presence of the superconducting leads (the "proximity effect" ) . 



A. Multichannel Luttinger liquid 



It has long been appreciated theoretically that the interacting Luttinger liquid provides a one-dimensional metallic 
counterpart to the higher-dimensional Fermi-liquid state. However, there are no fermionic quasi-particles in a Luttinger 
liquid and the elementary excitations are rather bosonic collective charge and spin fluctuations dispersing with different 
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velocities. One needs to take interactions into account right from the beginning by using the bosonization technique. 
The power of this technique resides in the fact that it allows to express an interacting fermionic Hamiltonian in 
terms of a noninteracting bosonic one. Correlation functions of fermion operators become correlation functions of 
(combinations of) bosonic vertex operators, which can be computed exactly. 

The bosonization procedure has been extensively reviewed in the literature ^ . Here we will only discuss in 
detail the particularities due to finite size and finite temperature, and how they depend on the boundary conditions. 
We start by noting that the fermion operators can be written in terms of bosons as 



La 

vr,J ^ i\/27[sinh 4,i , +cosh ^ „+a-[sinh tftj , +cosh „]] 



(32) 



where the ?7's are Klein factors obeying a diagonal Clifford algebra, ensuring anticommutation of different fermionic 
species. The index j = 1, ...N labels the channel, and cr = ± labels the spin. The interaction parameters t,^. ^ vanish 
in the noninteracting case, and are related to other conventional notations by 



= Kl 



(33) 



Repulsive interactions correspond to taking f < (K < 1), whereas attractive interactions correspond to taking f > 
{K > 1). In the following, we shall make use of both of these notations interchangeably. 

The Hamiltonian of the bulk interacting theory is that of free chiral left- and right-moving bosons for the charge 
and spin sectors: 



dxvl 



(34) 



The canonical commutation relations in the interval ]0, R[ are 



c^i^{x, t),d,<p{^{x', t)] = '-5^^'5ab5{x - x'), 



-S"'6atS{x-x'). 



(35) 



When working in the Hamiltonian formalism, our first objective is to provide adequate mode expansions for our 
fields. In a finite-size setup such as the one we are considering, the boundary conditions influence the particulars 
of this mode expansion: boundary conditions are in fact constraints on the fields that have to be properly taken 
into account. Unfortunately, in the interacting case, we cannot simply write down some interpolating boundary 
conditions like we have done for the fermionic fields in the noninteracting case. Instead, one has to impose some 
tractable boundary conditions, and use a perturbative setup around them. Two cases are detailed below, which we 
have already introduced in the noninteracting case: normal (eq. [t]) and Andreev (eq. boundary conditions. 



1. Normal boundary conditions 



The philosophy that we adopt for now is to start by imposing normal boundary conditions on the fermions (0) . In 
terms of the charge and spin bosons, this means that we fix (being careful with the zero modes) 



t) = ^{0, t) + 4>iR{0, t) = + n\] 



with n,m Cz 'Zi. The mode expansions satisfying these boundary conditions can be explicitly written down as 

<^iLi^,t)^^+Ili^^+C,ix,t), 



(36) 



(37) 
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with dynamical parts 



and commutation rules 



(38) 



(39) 



with the left and right zero modes commuting. 

The eigenvalues of the zero modes are set by the boundary conditions to relative values given by the topological 
number combinations (pq): 



Kr- 



2ai, Hi = (3i- al 

In terms of operators, the charge and spin Hamiltonians become 

HI - }^fv^ + ^\^ b b -^n^2^V6 b 

n^O n>0 



24i? 



(40) 



(41) 



Correlation functions can be obtained exactly in the finite-size geometry we are considering, even in the presence of 
a finite temperature. Our particular choice of fermionic boundary conditions yields a simple variation on the so-called 
open boundary conditions used in [ ^Of , where the computation of the two-point correlator is made. In the Appendix 
|A|, we have extended this result to the general case of the 2n-point correlation function. One thing to keep in mind 
is that we have computed canonical correlators, not (imaginary) time-ordered ones. The latter, however, can be 
obtained as combinations of the correlators we have derived. 

The first thing that we are interested in computing is the correction to the partition function coming from the 
boundary couplings representing the effect of the superconductors. We perform this computation starting once again 
from a real-time formalism. There are some subtleties when deriving a Matsubara-like formalism in the particular 
finite-size geometry that we are using: we devote Appendix ^ to this issue. 

The perturbing Hamiltonian, after making use of the boundary conditions, becomes 



Hi = H 



IL 



HlR, 



N 



N 



H,R = 2 2^ A'^ 



g»x+'*(fe^+'c^p)fl^'4^(^)^Jt| (i?) + h. 



(42) 



The correction to the partition function comes from the usual Matsubara term. The first order vanishes; after a bit 
of algebra, we are left with the second order ^-dependent contribution 



< p^^l 16C0SX ^l^R^ / d^i / 
ij = l -^0 -^0 



(1 - 5,,,)Gf\^i{Ti - T2),R)Gf\^t{n - T2),i?)+ 



-l-J,,,Gf'(-i(Ti-r2),i?) 



(43) 



where 3? denotes the real part, and the correlation functions can be computed using the formulas in Appendix 
After a bit of algebra, and some simplifications, we obtain for the two-point contributions the expressions 



a—c,s 

where the zero-mode part is 

,^£«lll^i(T^[fi + fcF^]kc)^4(r]f^|rl)+04(r^[§ + fc>i?]|T^)ei(rlf^|T]) 
Qi{-iT,R) = -e 3 



(44) 



02{Tj:k^pR\Ti)9,{O\Tl) + 03(T^fc>i?|T^)02(O|rj) 



(45) 
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The reader unfamiliar with ^-functions is referred to Appendix B. In order to make the bulk of our paper less clogged 
up, we have defined all functions and parameters in Appendix A. SufRce it to say here that the dynamical functions 
F are given by 

F (z) ^^(SK) (46) 



in which a is an implicit cutoff, i.e. z z ~ ia. The periods of the theta functions are 

(47) 



The phases are Pa/i = « . For the four-point function, we get 



G^)(-^r,i^) = A/^igf)/^, [Fl{-^v:T + R)]-''~^''^ (48) 



where the zero-mode part is 



-(4), . . ^ ^ Url[^ + fc>fl]|T^)^3(0|rl) - e^{Tl[^ + iepR]\Tl)e2{Q\Tl) 

A few comments might be in order here. First of all, the definitions of the functions used here are to be found in 
Appendix A. For a quick match with the usual correlation function behaviours, however, one just has to remember that 
the functions Fa{z) behave like cx z for zero temperature and infinite system size. Therefore, the scaling exponent of 
correlation functions is just (minus) the exponent of F in the expression encountered. This gives the usual bulk power- 
law correlation functions so typical in this kind of business. In particular, at the level of the dynamical contributions 
through the F functions, the charge and spin sector arc completely separated, a phenomenon which is usually called 
spin-charge separation. The two types of excitations, spin and charge modes, move at different velocities, seemingly 
splitting the electron in different parts. Notice the crucial fact, however, that the zero-mode part above does not 
separate at all in distinct charge and spin sectors. Via the zero modes, these two sectors remain completely entangled, 
unless we neglect all finite-size effects. In the limit of an infinite size, the zero-mode parts simplify to trivial phases. 
We are, however, interested in finite-size effects, so we keep these contributions throughout our treatments. This will 
turn out to be true of all correlators we shall compute in all cases: the dynamical sector readily factorizes into charge 
and spin sectors, but the zero-modes entangle the two in a complicated manner. 

The expression for the Josephson current to this order in perturbation theory is easily obtained by applying the 
derivative trick: 

lix) = 32sinx V A^A)?5i^ dr, H dr, [(1 - 5,,,)Gf ^ (-z(ti - ra), i?)Gf (-»(n - r2),i?)+ 



+(5,,,Gf^(-^(Tl-T2),i?) . (50) 



One thus obtains a sin^ behaviour to this order in perturbation theory, with a normalization that is bilinear in the 
boundary pairings, and contains a rather nontrivial function of the system size, temperature, velocities and interaction 
parameters. This reproduces the expected form of the current near the normal fixed point. Higher-order terms in 
perturbation theory could now be taken into account, which would start showing deviations from the pure sinx 
behaviour. Moreover, each term in the perturbation theory would be composed of some nontrivial function of x 
multiplied by a complicated integral over correlation functions, in other words a rather daunting function of system 
size. These would in turn conspire to make the current some scaling function of appropriate powers of the boundary 
pairings and the system size. We have not included those computations here, although they can in principle be 
performed using our formulas, given enough patience. 

One of the features of equation (^^ is that the current contains two different types of scaling functions of system 
size, with different exponents associated to them. In the limit of large system size and zero temperature, we can 
extract the leading behaviours of these, and write symbolically 
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The first type of scaling comes from Cooper pair tunnelling purely into one of the channels, whereas the second 
type comes from Cooper pair tunnelling into two different channels {i.e., the Cooper pair splits up in two single- 
electron components penetrating channels i and j). The particular form of these power laws depends crucially on the 
assumption that the system sits close to normal boundary conditions: if it is close to Andreev boundary conditions, 
the scaling is different, as will be seen later. One immediate thing to notice from these scaling functions is the expected 
fact that the current gets suppressed by increasing the repulsiveness of the interactions (that is, making smaller), 
whereas throwing in more attractive interactions enhances the current. This is very natural: we are dealing with a 
supercurrent, produced by pairing, and pairing works hand in hand with attractiveness. 
Another object of interest is the pairing operator, defined as 

= *i(a;,t)*f(x,t) (52) 

within the bulk of the system. The fact that the expectation value of this operator can be nonvanishing in the bulk of 
the system represents the effects of the superconducting boundaries: these produce the well-known "proximity effect", 
of which this is one manifestation. Using the normal boundary conditions, we can rewrite this pairing operator purely 
in terms of left-movers, and compute its expectation value to first order in perturbation theory. The two boundaries 
are independent to this order, so we take into account only the left one for simplicity. The result is 

<P«,.,0),.4A«£'.4,l-,,,)«{cf',-..,.)}»{cf(-„..)}-,,[«{G»,-.,.)}+C»,-^ 

The two-point functions in cither channel are here given by slightly more complicated functions (derived again from 
the formulas in Appendix A): 

Gl'\-^T,x)=N■,-'Qd-^r,x) [i^:(-z<r ~ x)] [F^(-^«lr -f x)] V:(2x)r ^ ^"^^f" (54) 

a—c,s 

where \F\ represents the modulus of _F, and the zero mode parts are again given in terms of ^-functions as 
^ _ . Wi 02{<[% + k^pR] + ^\rl)dM% + WVl) + UrlV% + fej^i?] + g|r-)g2(r]fg + i||r]) 

The four-point functions are of two types. This comes from the fact that, when writing the full pairing operator in 
terms of left-movers only, one obtains terms like 5'L|(a;)^'L|(a;) and 5'l|(— x)5'Lx(a;). These turn out to have slightly 
different x-dependencies. We label the first type -|-, and the second — . 
The first type of correlator is 

G%{-iT,x) ^ M:,^Q%l{-iT,R) \F,{~iv,T ~ x)\-^^^'~'''^ \F,{-iv,T + xf~'~''' |Fe(2a;)r^'"h2«. (gg) 

. 02(r^[^ + fc^i?] + ^|r^)e3(0|T]) + 03(r^[^ + -f ^|T^)^^2(0|r]) 

,x) — e 

The second type is 



with 



' ' 02(r^A:],i?|T^)03(O|Tj) + 03(r^fc>i?|T^)^2(O|rJ) ' ^ ^ 



G^^l{^tT,x) - gL'i(-*T,i?)|i^,(-it;cT + a;)r2-"'n^c(2a;)|4^"'nF42a;)r^'=''^ (58) 

with 

^(4), . , ^2(t^[^ + fc|^i?]|r-)g3(5f kl) + g3(T^[f + fc>i?]|r^)^2(5f |r]) 

^-'"^ ' ' ^?2(r^fc>i?|T^)03(O|r|) + 03(T^A:^i?|r^)e2(O|Tj) ' ^ ^ 



One therefore sees that the three terms in (^3|) have different scaling behaviour when interactions are added to the 
bulk of the system. In particular, the behaviour of the pairing expectation value depends on whether we are looking 
at channel-diagonal or cross-channel pairing. The results can be summarized by the scaling laws 

In these, x represents the distance perpendicularly away from the superconducting boundary. The cross-channel term 
has a decay law that depends on the interaction parameters in both channels, whereas the channel diagonal term has 
two different behaviours, with different exponents depending on the nature of the operator involved. These results 
are collected in simplified form and discussed further in section fv\ . 
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2. Andreev boundary conditions 



The other case that we consider, which leads to tractable calculations, requires extremely good coherence between 
the original bulk superconductors and the one-dimensional channel. This does not mean a perfect contact in the usual 
sense of the word, however. We invite the reader to consult reference |2^, where basic considerations on the level of 
a simple lattice model point to the fact that there exists an optimization phenomenon, by which it should be possible 
to achieve perfect Andreev reflection at the boundaries by tuning available parameters. 

We therefore want to reproduce the above results in this case. Since the specifics of the boundary conditions are 
crucial to the physics of the system, physical quantities like the Josephson current and the superconducting order 
parameter profile will be different from the previous results computed around the normal fixed point. We therefore 
start by imposing Andreev boundary conditions on the fermions, eq. (14). The right-moving fields can then be 
described in terms of left-moving ones via 



with periodicity conditions 

^La{x + 2R) = -e'^^Laix). 

In terms of the charge and spin bosons, this means that we fix 



(61) 
(62) 



TT 



0i(O,t) 



0i(i?,i) 



2K 



2K 



7 

3 ^ T 



-1 = 0^ 



(63) 



with n^m Cz Z. The mode expansions satisfying these boundary conditions are the same as for normal boundary 
conditions for the spin channel. For the charge channel we have instead 



with dynamical parts 



m- 



2R 



(64) 



(65) 



ne2 



and commutation rules 



id'''. 



[6; 



cn' cmi 



nS, 



(66) 



with the left and right zero modes, and charge and spin modes commuting. 

The eigenvalues of the zero modes are the same as before for the spin modes, and for the charge sector they become 



ui = Pi 



(67) 



The fundamental difference between normal and Andreev boundary conditions is that for the Andreev case, there 
exists a Josephson current and a nonzero pairing operator to zeroth order in perturbation theory. That is, we can 
first of all compute the partition function in the Andreev case, which carries a phase dependence giving nonzero I{x) 
by using the derivative trick. 

The partition function is given hy Z — Tr[e~^^ j^'^'"'^^''^] where is the charge or spin Hamiltonian for a given 
channel ( eq. ^l|). Since the Andreev boundary conditions as well as the interactions are channel diagonal, the total 
partition function is a product of the partition function for each channel. Therefore the total current will be the sum 
of the contributions of each channel. In other words, we need to compute simply the partition function for one LL 
coupled to Andreev boundary conditions. This computation was already done by Maslov et al. in reference [|l^. Here 
we only sketch the derivation. Since the Hamiltonian can be written in terms of a zero mode part plus a dynamical 
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part that commute with each other, the partition function becomes a product Z = Zzm('x) x Z^yn where only the 
zero modes part depends on x- Using the exphcit expressions for 11^ and given in eqs. and ( |67| ) respectively, 
and following the same steps described in Appendix A for the computation of the correlation functions, we obtain 

Z.rn = e'^[Mrc^\rc)MO\rs) + Mrc^\rc)MO\rs)] (68) 

where Tc — if3^^^^ — . We have dropped the channel index for simplicity. The current can be readily obtained by taking 
the derivative with respect to x of the above expression for each channel, and adding over the channels. Since the 
general form of the current for arbitrary temperature and length, involves derivatives of the i?2,3 functions, it can not 
be written in a closed form. The asymptotic values of the current for a given channel are: 

I'ix')=e^^ (69) 



R 



at low temperatures (^^ >> 1), and 



P{x^)^SeTe '\ "W sin(x^) (70) 

at high temperature << !)• This expression reflects the fact that at high temperatures, repulsive interactions 
{Ki — e?^' < 1) suppress the Josephson current [^s). 

Here as well, the channel diagonal component of the pairing operator (eq. (p2[)) can be computed exactly , without 
the need for perturbation theory. The general expression for P"{x,t) can be obtained from the equations in the 
Appendix A by taking the coordinates of both fermions at the same point. Here we only write the asymptotic values 
of the order paramater at high and low temperatures. 

At low temperatures we obtain 

<P-(x)>«e^t-(sin(^))"^^""^^"^ (71) 

In the limit of infinite length this correlator behaves like < P"{x) >~ with 7 = {e^^^" + e'^^=)/2. This result was 
already derived by Maslov et al. in reference [Q. 
At high temperatures, the order parameter becomes 

< P"(a;) >« e^'^f (^^+'^^ (72) 

We emphasize again that the two results above are completely nonperturbative, since the partition function and 
correlators can be calculated exactly at finite temperature in this finite-size geometry for what we have termed Andreev 
boundary conditions (eq.(p^)). It is however possible that in the physical system described by the multichannel 
Luttinger liquid , the coupling to the superconductings leads is represented by non diagonal matrices A/j ^ . This is 
certainly the case if the distance between the channels in the Luttinger liquid is smaller than the coherence length 
in the superconducting lead. In this case, it is necessary to take into account the off diagonal contributions of the 
boundary action, ie, those terms corresponding to l 7^ for * 7^ J- This can only be done perturbatively, in an 
analogous way to the one described in the previous subsection. We do not include explicity the calculations here, 
but the contributions of these perturbations to the scaling behaviour of the superconducting current and the pairing 



operator are summarized in section IV. 



Carbon nanotubes 



Carbon nanotubes (CN), long, thin cylinders of graphite discovered in 1991 by S. lijima jsij, are large molecules 
that are unique for their size, shape, and remarkable physical properties. In fact, they can be metallic depending 
on the wrapping of the tube (|, ||, 32 . In the following we will be interested in those CN which are metallic and 



are known as "Single walled nanotubes" (SWNT). They exhibit unique quantum wire properties that derive from 
the tubes' nanometer scale diameters in combination with the special electronic structure of graphite. Because of 
the quantization of circumferential modes, the tube's electronic states do not form one wide electronic energy band 
but instead split into one-dimensional subbands with band onsets at different energies. These subbands are widely 
separated in energy, on the scale of 1 eV (much larger than the room-temperature thermal energy fc^T, about 0.025 
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eV). Only two of the one-dimensional subbands cross the Fermi energy in metallic carbon nanotubes; all the current 
through such tubes is therefore predicted to be carried by only this pair of subbands. Thus, a single wall carbon 
nanotube can be described by a pair of Luttinger liquids whose velocities are equal to each other, and that interact 
through Coulomb interaction. We shall therefore consider only two conducting spinful channels, i.e. (taking spin and 
chirality into account) eight bosonic fields with a — L, R, a = ± and s = 1,2. The main difference between 
the effective description of a carbon nanotube and the previous case of multichannel LL, is that for CN the bulk 
interactions are not channel diagonal. This is due to the fact that the Coulomb interaction couples only to the total 
charge density ^ while the spin sector remains unaffected by the interactions. Therefore, we move to a basis of 
total and relative bosons: 

•P;.a^^i<t>la+r<t>la)- (73) 

with T — ±, fj, = c, s, and a — L, R. In terms of these operators, the Hamiltonian becomes (after taking the interactions 
into account nonperturbatively) that of a set of free bosons, 

H=Y.ll f '^^"A^ [id-^lL? + {d..4'lR?] (74) 

^=C^S T = ± 

which we can now use to represent our fermions with the bosonization formulas 

La y^'lLa^ 

In the above, is the interaction parameter representing the Coulomb interaction that involves the total charge of 
the system (in both channels of the nanotube). 



1. Normal boundary conditions 

As was the case with the Luttinger liquid, in order to set up our computations correctly in the finite size geometry 
that we are considering, we have to specify tractable boundary conditions. We start with normal boundary conditions 
for the fermions, which correspond in the bosonic language to taking 











± 




T 





^^^^^ ^ + ] ^ «^ = ^1 

<l^t{R) = + "^f + 1 ± 1 - ^] ^ Pc. ^f{R) ^^H- "^f ] ^ /3?- (76) 

where — Kc, K~ = 1, = ± with independent integers, and kp — kp ± kp. The mode expansions for 
charge and spin operators coincide with the ones given in eq. (|37|). 

The second-order correction to the partition function can be calculated in a similar manner as for the multichannel 
Luttinger liquid. After a bit of algebra, and again making use of the correlators in the Appendix A, we find the 
X-dependent contribution 

(pf^)x = 16cosx J2 ^i^R r^^i r'rfT25R[(l-%)Gi^^(-«(ri--T2),i?)+5„G{|"^(-z(Ti-T2),i?)l , (77) 

^,j = l,2 Jo Jo ^ J 

with correlators 

g[|^^ {-it, R) = P+Q(ui)D|| (-ir, R), 

Gf{-iT,R) = P+Q{n\\)D^{-iT,R), (78) 

where we have defined P+ = e'^I^^*^ '^"^ (we remind the reader that Kc — e^^" is the only nonvanishing ^ for the 
nanotubes). The dynamical parts are for both correlators equal to 

D\\{-iT,R) = {F,{-iTv,^R)\-^'^'^ {F,{-iTv, + R)\-^ , 

D^{-iT, R) = [F,i-iTv, + RT"'''^ [F,{-iTVs + R)r' , (79) 
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independently of the channel index, and the zero-mode part is given in terms of the four-dimensional 0-function 
(defined in the second appendix) 



g(u) = e 



_ i{2kFR-TT)T^T/P 



9{u + ^{2kFR ~ Tr)m\n) 
6'(^(2fcFi?- 7r)m|rj) 



(80) 



with vectors 



/ iT[e 



TT 



\ iT[e 



2R\ 



2R J 



IT 

AR 



iTVc[e 
iTVc[e 



1] 
1] 



V 



iTVc[e - IJ 
iTVc[e~^^'' — 1] 



2R \ 
2R 

J 



(81) 



and m = (1,1,1,1). The interaction-dependent period matrix Jl is defined in equation ( A27 ). We also remind the 



reader that Tc 



^ ^""fl • '^'^ obtain the above results we have used that G,, ^ 



(4)12 _ ^(4)21 



Gil ^ due to the symmetry 



properties of the matrix Jl. 

For the Josephson current, this then implies that, to second-order in perturbation theory, we find the expected sinx 
behaviour. More precisely, we get 



I{x)^F(P,R)sinx, 



(82) 



with a normalization amplitude function that depends on all the parameters of the system, and is given in terms of 
the integral expression 



F(/3,i?) = 32^A^AHi 

P Jo 



(1 - 5,,)G'-l\-iiri - T2),R) + ,5yG[f^(-i(ri - ra), i?) 



(4)/ 



(83) 



As was the case for the multichannel Luttinger liquid, the scaling behaviour of the current splits up into two different 
types, depending on the way the Cooper pair penetrates into the conducting channel. Here, for the case of normal 
boundary conditions, the scaling behaviour is the same for channel-diagonal and cross-channel pair penetration (i.e., 
doesn't depend on whether the pair penetrates into one channel only, or splits between the two): 



/(x)^(/||+/±)i?- 



(84) 



The physical interpretation of this is easy: as Coulomb interactions involve the total charge density in both channels 
together, it doesn't matter how the pair tunnels in, as the interaction cost is the same in both cases. 

For the pairing operator defined in equation (52), a first-order perturbative computation similar to that in the 
Luttinger liquid case yields 



(P^^(x,0)) 









(!-%)( 


Jo 



(1 - %) (G^ll{-iT,x) + G^^l(-ir,x)) + (G\^li~iT,x) + Gf}i-iT,x))] , (85) 



with correlators splitting up into two different classes (channel-diagonal or cross-channel) and again two subclasses 
(±, depending on the particular operator combination; we refer the reader to the note at the end of the computation 
of the pairing for the multichannel Luttinger liquid with normal boundary conditions): 



.- (l±l)7rx 



G'^ll{-iT,x) = Q{uj_±)e'^' ^r"" Fj_±{-iT,x), 



Gl^h-iT,x)=Qiu«±)e'^^F«±i-iT,x) 



(86) 



with dynamical parts given explicitly in terms of powers of _F-functions (note in particular that if charge and spin 
velocity coincide, then the subsets ± become identical) 

\Fc{2x)\-'''^ 



F±+i- 


iT, 


x) 


= [Fc{- 


irvc 


-x)] 






[Fci- 




+ x)Y^^'' 


F^-i~ 


IT, 


X) 


= \Fc{- 


irvc 


+ ^)l 




iFsi^ 


-ITVs 


+ ^)l 


-'\Fc{2x)\ 




IT, 


X) 


= [Fc{- 


irvc 


-x)] 






[Fci- 




+ x)Y^^-' 




IT, 


X) 


= \Fc{- 


irVc 


+ ^)l 




iF^i- 


-ITVc 




-'\Fc{2x)\- 



[Fs 



\Fc{2x)\- 



x)y 



[Fc{-ii 



(87) 
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The zero-mode parts are here given by equation (g^), and the distance vectors appearing in the multivariable 0- 
functions representing the contributions of the zero modes are defined by 



/ 



U_L± 



TT 



vs] ±2x\ 



iT\e 



U||± 



TT 



2x j 



«TWc[l 

irwn[l 



± 2a; \ 



2x 



zTWc[— 1 + e ^^^] 



y iTUc[— 1 + e 



We therefore again find three different scaling behaviours for the pairing operator, as in the multichannel Luttinger 
liquid case. The scaling of the pairing operator expectation value again splits up in two different types of behaviour, 
which we can summarize as 



(89) 



First of all, we note again that the scaling of the + operators coincide, as does the scaling of the — operators. That 
is, it's not important (as far as scaling is concerned) how the Cooper pair penetrates the nanotube. These results are 
summarized and discussed in section IV. 



2. Andreev boundary conditions 



For the relative bosons suitable to describe carbon nanotubcs (73), Andreev boundary conditions can be written as 



4>fiR) = ^^-^K + "^f - 1 T 1 - ^] ^ 0t (i?) = ^[mf - mf + i ± i] ^ (3f, (90) 

where = n^. ± with independent integers, and = ^ X^- The mode expansion for charge and spin 
operators coincide with the one given in eqs. (^) and ( |3^ ) respectively. 

In order to compute the superconducting current, we consider the partition function. Once again, the only part 
of the partition function which is x dependent is the zero mode part Zzm- Since the two channels constituting the 
nanotube are coupled to the superconductor leads at the boundary in the same way, the two phases must satisfy 
X^ = x'^ ~ X, E^J^d the current is determined by Zzm{x) 

Z^Ux) = e'^0(^m|Jl) (91) 

where fc = iP^^^ — , and rii — (1, 1, 1, 1). The multivariable theta function 9 and the matrix period $7, defined in 
Appendix B and A respectively, appear naturally when computing correlation functions for carbon nanotubes, where 
the electronic interactions are not channel diagonal. 

The total current is obtained by deriving In Z^m with respect to the phase x- Its asymptotic values are in this case 
given by the expressions 

/(X) = e^^ (92) 

H IT 

at low temperature >> 1), and 



I{x) = leeTe-'^t^" sinx (93) 

at high temperature (^ << 1). 

A physical realization of the coupling between a nanotube and a superconductor would however produce boundary 
pairing amplitudes A*-' which should be approximately equal for all the channels (that is, even for the off-diagonal 
components). This is a fundamental difference between the present case and what one could expect from a two- 
channel Luttinger liquid: in the latter case, the off-diagonal pairing components could be taken to vanish if the 
contacts between the original superconductor and the two Luttinger liquids were separated by a finite distance. The 
diagonal components of the boundary pairing are easy to treat nonperturbatively as seen above. In view of the last 
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comments, however, we also want to include another contribution to the partition function coming from the channel 
off-diagonal boundary pairings. The perturbing Hamiltonian is given by 



Hi = H 



IL 



H 



IR, 



(94) 



where 



HiL = 



1 



Al2,T,lt^n^,T,2t 



(0)*f (0) + Af^f (0)^1^0) 



Hu 



1 r 
2 



h.c, 
- h.c. 



(95) 



The first non vanishing correction is the second order one, which is given by 



< Pelf >x~ 



2AlAr^ dTi I dre 



i2kF R—ifc^T 



-IVcT 



R)]-"^ [F,{~iVcT - R)] [F,(- 



where 3? denotes the real part, the vectors U| and uj^ are defined by 



U| = 



4i? 



V 



~i{-Vre'^'^' + Vs)t ^ 
-2i? + i(wce^^'= +Vs)t 
-2R + ilvcC^^" +Vs)t 

-i(-Uce^^- + Ws)t / 



TT 



iVsT + R)] ' 



/ -2i? + i(uce2«<= +Vs)t 
-ii-Vce'^^" + Vs)t 
-i{-Vce^^'' + Vs)t 
\ -2R + iivcc'^^-^ +v,)t J 



(96) 



(97) 



respectively, and the correlation functions were computed using the formulas in Appendix W. The correction to the 

(2) — 

superconducting current will be given by /'■^-'(x) — ~^~^^^§^^^- Even though wc can not give a closed expression 
for this correction, we can study some simple limits. In particular, we can show that in the limit of low temperatures 
/(^^ (x) "^ill depend on the length of the system as times a complicated function of x, representing deviations 
from the sawtooth form of the perfect current. Thus, the off diagonal contribution to the superconducting current 
introduces a dependence on the length of the system which is interaction dependent. Due to the general structure 
of the correlation functions given in Appendix ^ , it is simple to see that higher order terms in the perturbative 
expansion will contribute with higher order powers of to the superconducting current at zero temperature. 

The order parameter can be computed exactly if we only consider channel diagonal Andreev boundary conditions. 
In this case, using the boundary conditions (eq. |6^) we can express the pairing operator (eq. (p3)) as 



< P"(x) >= 



< 



> 



C-frh\n) 



[Fc{-2x)]- 



m2x)Y 



[Fsi-2x) 



(98) 



where = ||(1 + a,l - fi, 0, 0) and = ||(0, 0,1 -f ct, 1 - a). The contributions to < P"(a;) > of the form 
< ^^|(x, t)^'^|(a;, i) > as well as the off-diagonal components of the pairing operator (< P^'^{x) >,< P'^^{x) >) 
vanish unless we take into account off diagonal terms in the boundary couplings. 
At low temperatures we obtain 



< P(x)" >w e'^x(£ 



R 



(99) 



In the limit of infinite length this correlator behaves like < P{xy 
At high temperatures, the order parameter becomes 



with 7|| = (e-2«c +3)/4 



< Pixy >w e' 



(100) 



However, this is not the physical case since the LL constituting the CN couple to the boundary in the same way, 
introducing Andreev reflections between the two channels. Thus, we need to include perturbative corrections around 
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the channel diagonal case. The first order correction due to the off-diagonal Andreev terms (eqs. (|95[)) vanish for the 
diagonal components of the pairing operator < P" >, but give a finite contribution for < P^^ >, and < P^^ >. Since 
to this order the two boundaries remain uncoupled, for simplicity we only give the correction due to the left boundary 



with correlators 

where the distance vectors are defined by 



/3 



dT 



(101) 



(Ifmin) 



[Fc{iTVc - x)Y 



\Fci2x)\' 



[Fs{ 



IT 

AR 



~tT[e--'=Vc - Vs\ 
-iT[e^^"Wc + Vs] + 2x 
-irfe^^^Wc + Vs] + 2x 

V. 



irle^^^Vc - Vs 
irle'^^^Vr + Vs] 4 



2x \ 



2x J 



(102) 



(103) 



Thus, there is only one type of non vanishing contributions to this operator (to first order in perturbation theory), 

namely, what we have called G^'. The other terms present in the pairing operator of the form ('I'^|(— x, t)^^'| (2;, i)) 
do not contribute to the expectation value of the off diagonal pairing operator to this order, due to the particular 
form of the boundary perturbation. 

In general, this expectation value will be a complicated function of the distance to the superconductor x. However, 
in the limit of zero temperature, and for very large systems, we can show that (P^^(a;, O))'"'^-' ~ with jj_ — 
(e-2«= +3e2«=)/4. 



IV. SUMMARY AND DISCUSSION 



Up to this point, we have presented most computations without making any approximations with respect to system 
size or temperature. This results in rather complex expressions for correlators and their related physical quantities 
which are challenging to interpret straightforwardly. We have done this because it is feasible, and represents a very 
good starting point for seeking a more precise exact solution to the type of problem that we have considered. For 
example, the exact correlators that we have computed would have to be reproduced by any other method claiming 
to be capable of dealing with similar systems (we are thinking here of for example massless limits of form factors for 
the sine-Gordon model). 

In order to make the important results more accessible, and to paint a more transparent picture of the systems we 
considered, we provide here a short summary of our results. We remind the reader that the two systems that we have 
dealt with, namely multichannel Luttinger liquids and carbon nanotubes, have different bulk interaction structures. 
We have considered charge and spin interactions which are channel-diagonal for the Luttinger liquids, whereas the 
interaction for the nanotubes is of the form of a Coulomb term involving the total charge of the system. This has an 
influence on the specifics of the bosonization procedure, that was discussed in the bulk of our paper. 

We have identified two different fixed points for our theory: one which we call the normal fixed point, and the other 
the Andreev fixed point. The normal fixed point is associated to the normal boundary conditions (^, and similarly 
for the Andreev fixed point (p^. At each of these fixed points, it was shown that it is possible to do the bosonization 
exactly in finite size. All correlation functions are then readily computable, and we listed all correlators for all cases 
in Appendix A. 

We have presented explicit computations for physical quantities like the Josephson current, and the pairing order 
parameter. We summarize here these results, and perform a stability analysis for the fixed points we have identified. 

First, we define /y and I± as the contributions to the Josephson current from channel diagonal and channel 
off-diagonal correlators respectively. Then we define the set of operators whose expectations values were already 
computed, and whose scaling needs to be assessed. The channel-diagonal and channel non-diagonal pairing operators 
are given by 

P||+ = P||_=r^^(x)*^,^(a:), 

Px+ = ni{^)n^i^)^ p^- = nii^)K^i^)' i^^j) (104) 
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both of which occur in the physical pairing operator involving the full fermion. For completeness , we also define the 
back scattering bilinears 

= ^%{x)^i,{x), 14_=*t(^)*L(^). (105) 

In Tables | and || we summarize the results obtained for the Josepshon current and pairing operator in the limits of 
zero temperature and infinite size. For the diagonal and off-diagonal components of the Josephson current we list the 
exponent that determine their dependence on the inverse of the system size, . We also include the decay exponent of 
the pairing operators as a function of distance away from the boundary. We list the results obtained considering the 
first non-vanishing contribution of the expectation value of any given operator, using the appropriate perturbation if 
necessary. 



TABLE I; Scaling exponents (Luttinger liquids). 



Operator Normal Boundary Conditions Andrcev Boundary Conditions 

1\\ 2e-2C _ 1 i 

P||,+ i[e2C + 3e-2«c] _ 1 1 [3e2C + g-^C] _ i 



TABLE IL Scaling exponents (nanotube) 



Operator 



Normal Boundary Conditions 



Andreev Boundary Conditions 



+ 



/|| 

I± 

P|l 
P||.^ 

P±.- 



e-2«c 



3e-2«c] 
+ 3e-2«. 



1 



g2Ce 



2Cc , 

i[3e2«=- 



e-2«. 

-3] 
e-2C. 

-3] 



As it is clear from the previous section, in the case of Andreev boundary conditions, the results listed for /|| and P||,- 
are exact, I± was obtained to second order in perturbation theory, and P||,+ as well as P_l,± were obtained to first order 
in perturbation theory. In all these cases the perturbative Hamiltonian corresponds to Andreev backscattering only. 
For carbon nanotubes and Andreev boundary conditions, we included in the bulk of the paper only the off-diagonal 
terms of the boundary perturbation. In this case < P||,+ >= 0. However the result listed in table || shows that this 
operator acquires a finite expectation value if the boundary perturbation includes diagonal terms. Finally, it is simple 
to see that the second order contribution to Pj_,- vanishes due the particular form of the boundary perturbation. 
Therefore, althought we include its naive exponent here, it should be remembered that in our setup, this vanishes. 
The scaling would then be determined by the first descendant of this operator. 

These tables are simply a convenient way of viewing the various dependencies on the interaction parameters. The full 
functional expression for the operators' expectation values, however, is much more complicated, containing nontrivial 
functions of the temperature, size of the system, velocities and interaction parameters. The results are too bulky to 
summarize here more precisely than the table above, and we refer the reader to the previous sections for the exact 
formulae. 

What we have to do next is c heck t he s tability of our fixed points by computing the effective dimension of the sets 
of operators given by eqs. (104) and (105), when bulk interactions are introduced. We provide four tables giving the 
scaling dimensions of the pairing and backscattering operators around both the normal and Andreev fixed points for 
the (multichannel) Luttinger liquids and nanotube. These exponents are all simplified data coming from our exact 
correlators, and they are valid for zero temperature in a large system. In the Luttinger liquid case, the interaction 
parameters carry a channel index, whereas we have considered here only one nanotube, with one charge interaction 
parameter. 



21 



Three exponents are given, corresponding to the various regimes of the correlators that can be encountered. A is 
the usual bulk scaling exponent, valid away from the boundaries. Aj^ is the scaling exponent close to one boundary, 
with argument separation running perpendicular to the boundary. Ay is the scaling exponent for correlators where 
operators both live close to the boundary, but with a finite time difference (the reader unfamiliar with these various 
exponents is invited to see e.g. [? ] for a clear explanation). 

To examine criteria of relevance in the RG sense in our theory, the exponents to look at are in the A|| columns. 
That is the case in our framework, since our boundary perturbations include only fields sitting on the boundaries, 
integrated over an infinite time interval; in principle, we should perform a complete RG analysis, including even 
terms which represent the appearance of superconducting order within the one-dimensional channel, away from the 
boundary. We don't to this here, as these contributions turn out to be much smaller. It is readily seen that, for 
normal boundary conditions, the pairing perturbations are irrelevant (Ay < 1) for repulsive interactions for both the 
multichannel Luttinger liquid and the carbon nanotubes. On the other hand, the diagonal backscattering term is 
marginal (the cross-channel backscattering is either irrelevant (Luttinger liquids) or marginal (carbon nanotubes)). 
Putting attractive interactions in the bulk, however, destabilizes the normal fixed point: in that case, the operator 
that drives the system away is the channel-diagonal pairing operator. One therefore expects, on general grounds, in 
the case of attractive bulk interactions (and supposing that our parameters are initially tuned to put our system near 
the normal fixed point), to witness a flow of our theory towards Andreev boundary conditions, corresponding to an 
optimization of channel-diagonal Andreev reflection at the boundary. 

At the Andreev fixed point, the operator scaling dimensions are different, and are listed in the last two tables. 
We stress this fact again: the particulars of the boundary conditions influence in a nontrivial way the scaling of the 
system near the boundaries, and near the Andreev fixed point, the dimensions of the operators differ from those to be 
found near the normal fixed point. The situation is now reversed as compared to the scaling around the normal fixed 
point: most operators are now irrelevant for attractive interactions, with a few of them being marginal. If we were to 
turn the interactions in the bulk to the repulsive regime, we would see the system flowing back to the normal fixed 
point. We therefore recover a similar situation to the one in |2^, in which the normal fixed point was the fixed point 
associated to repulsive bulk interactions, and the Andreev fixed point was associated to attractive bulk interactions. 
Our analysis here is rather short and simplistic: we could ask e.g. what other possibilities there might be for fiows in 
the system, in particular by playing around with the interactions in the spin channels. A treatment similar to that in 
would then be obtained. We leave this for future considerations. 
The perturbative setup that we have presented here lends itself ideally to further extensions and refinements. First 
of all, the formulas we have presented for the noninteracting case would be appropriate to investigate the effects of 
voltages simultaneously present with various backscattering potentials, or could easily be modified for the addition 
of ferromagnetic couplings. Other possibilities more directly related to what we have done in the interacting case 
would be, for example, to take into account perturbatively the effects of quasiparticle penetration within the original 
superconductors. We have neglected those according to the logic of and this should be correct at extremely 
low temperatures. However, finite temperatures make quasiparticle penetration possible if the energy scale starts 
approaching the (original, bulk) superconducting gap. This could again be done within our framework, by replacing 
the boundary action with some other action containing retarded effects. It is also important to realize that the class 
of problems that we have dealt with here forms a good basis from which to study other geometries. In view of the fact 
that experimental work is addressing these issues on a frequent basis, we think that there exist many opportunities 
for simple adaptations of our method. One could consider things like nanotube crossings of various kinds, tunnel 
junctions, etc. We leave such issues for later work. 



TABLE III: Operator scaling dimensions (Luttinger liquids with normal boundary conditions). 

Operator A A_l Ajj 

P||+ cosh(2e) i[e2«e +3e-2C] e-2«e 

P^^± i[cosh(2e) + cosh(2e) + j)] ^^[e^C + 3e-^il + e^?: + Se-^H + (^ ^ j)] l[e~^il + e'^^l + (z ^ j) 

1^1, _ i[e2«i+e2f:] i[e2C+e2«:+2] 1 

yj_,± i[cosh(2C) + cosh(2e:.) + (* ^ j)] j^[e2C + 3e-^C + + 3e-2«: + (z ^ j)] \[e-^C + e^^c: + (, _ 
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TABLE IV: Operator scaling dimensions (Luttinger liquids with Andreev boundary conditions). 

Operator A A_l Ajj 

cosh(2e^) i[3e2C+e-2C] e^C 

Pli- i[e-2«^+e2«:] i[e-2C+e2C +2] 1 

P^,± i[cosh(2a)+cosh(2e^) + (z^j)] 3L[3e2«:+e-2fc + e2€+3e-2f:+(i^j)] ^[e^Sl + e-^^l + {i ^ j)] 
V\\,+ 1 1 1 

yj_',± i[cosh(2^j) + cosh(2e) + (i-j)] ^[3e2C+e-2C+e2C+3e-2S^ + (i^i)] i[e2C + e'^?: + (i ^ j)] 



TABLE V: Operator scaling dimensions (nanotube with normal boundary conditions). 

Operator A Aj_ Ajj 

i[cosh(2^,) + l] i[e2fo+3e-2C.+4] ^[e-^^^ + l] 

P||,_ i[e-2«»+3] |[3e-2«=+5] ife-^^-^ + l] 

P±,+ i[cosh(2Cc) + l] i[e2«=+3e-2«c+4] ^[e-'«" + l] 

P±,- i[e-2«°+3] |[3e-2«c+5] 5[e-'«= + l] 

1 1 1 

Vii,^,- i[e2«=+3] |[e^«^+7] 1 



TABLE VL Operator scaling dimensions (nanotubes with Andreev boundary conditions). 



Operator 



Ai 



^11,- 
P±,+ 
P±,- 



i[cosh(2Cc) + l] 

i[e-2«= 



i[3-2«. 



3] 

i[c0sh(2^e)- 



+ 3] 



i[e2«c+3] 



7] 



|[3e2«=+e-2C= 
i[e-2«.+7] 



i[3e2«c+5] 



+ 4] 

+ 4] 



i[e2«c + l] 
1] 



2 
1 



APPENDIX A: CORRELATION FUNCTIONS 



The usefulness of bosonization really manifests itself in the computation of correlation functions. In this appendix, 
we compute multipoint correlation functions, in finite size and at finite temperature, for all the cases of interest in 
this work. 

The first requirement when dealing with a finite size is to specify the boundary conditions used in the computa- 
tion. We consider here normal and Andreev boundary conditions introduced in the bulk of the paper. We present 
computations both for the Luttinger liquid, and for the carbon nanotubes. Although two-point correlators for such 
systems are already known in finite size and at finite temperatures, and multipoint correlators in the bulk are also in 
the literature, we wish here to present formulas for any multipoint correlator. 



23 



1. Luttinger liquid, normal boundary conditions 

Let us consider the problem of computing correlation functions in finite size, at finite temperature, for normal 
boundary conditions. The general correlator that we want to compute can be taken to be 

n n 
i=l 3=1 

since all others can be obtained by simple manipulations (here, we mean in particular correlators with a different 
ordering of the operators, which are easy to recover from the formulas we present; see later for a better 

explanation). Correlators involving right movers can be recovered by rewriting the right movers as left movers 
using the appropriate (here, normal) boundary conditions. Here, the angular brackets denote an average using the 
Hamiltonian (|4l|), and the subindex TV denotes normal boundary conditions. 

For simplicity, we have suppressed the channel index. We remark that when dealing with multichannel LL, all the 
correlators needed to compute various physical quantities can be reduced to the above one for a given channel , or 
products of them evaluated in different channels. 

The first step in the computation is to substitute the bosonization formulas (|3^) for the fcrmions, and to make 
use of the normal boundary conditions, mode expansions and Hamiltonian (equations ( ^6| ) to (pl|)) to reexpress the 
correlator as a product of zero-mode and dynamical contributions, averaged over a free theory. We write this as 

^{l"r}({^' "}) = ^Y. '^.T, rJ{a^T} {{x, t] y, u}) K i^^-r}{{x , t] y, u}). (A2) 

Let us first consider the zero-mode part. Using the Campbell-Baker-Hausdorff formula to combine exponentials, 
we can reexpress this contribution as 

J{a,T}{{x,t;y,u}) = Pc{{x,t;y,u})Ps{„^r}{{x,t;y,u})Q{„^r}{{x,t;y,u}) (A3) 

with phase functions 

/ n n \ 

InPc = I n'^lvce^'^^-' {ti - m) +Xi-yi\- ^ [vce^'^^''{U - tj + Ui - uj) + Xi - xj + yi - yj] | , 



i=l i<j = l 



^'aPs{a,T}i{x,t]y,u}) = — — (^a;]^[a,{vse '^^'U + x,) - Ti{vse u^ + yt)]- 



\ I l—l 

\ 

- ^ [aiaj{vse {U - tj) + Xi - xj) + TiTj{vse '^^^ {u^ - Uj) + yt - yj)] (A4) 
i<i=i J 

and zero-mode expectation value 

Q{a,r}{{x,t;y,u}) = (e*v^mSr=l["oe-«'^(t,-«0-He«'^(x,-y,)]g»v^SSE^^Ji,,e-«Ma,t.-r.«.)+e«M<T,x,^ ^^g-j 

This average is readily computed in ter ms o f functions (see next appendix for definitions and useful formulas) . The 
final result is given below in equation (A16). Note the crucial fact that the spin and charge sectors remain entangled 
in this contribution. 

The dynamical part, on the other hand, readily factorizes in charge and spin sectors: 

K{cr,r}{{x,t;y,u}) = Kc{{x,t;y,u})Ks{a^^}{{x,t;y,u}). (A6) 

For clarity, we will here present details of the computation of the spin part only (the charge part doesn't depend on 
the cr, r subindices and is simple to recover from the formulas below). As an expectation value, this is 

n n 
1=1 1=1 

Since the average is over a free bosonic field, we can make use of the identity 

([]e^0=e<[^-<.^-^^ + 5i:.^?]> (A8) 
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to rewrite the above expectation value as 
lnK^[„^^y{{x,t;y,u}) = 

= 2n cosh^ ^a[^aiTjBs{xi,ti;yj,Uj) - ^(JiUjBs{xi,ti;Xj,tj) -^TiTjBs{yi,Ui;yj,Uj)] + 



+27rsinh^ S,s^[^(TiTjBs{-Xi,ti] -yj,Uj) - ^aiajBs{-Xi,ti; -Xj,tj) - ^nTjBs{-yi,Ui; -yj,Uj)] - 

i,j i<j i<j 

-7rsmh2^s['^aiTj{Bs{xi,ti;-yj,Uj) + Bs{-Xi,ti;yj,Uj)) - 
-'^(^i(^j{Bs{xi,ti; -Xj,tj) + Bs{-Xi,ti;xj,tj)) -^TiTj{Bs{yi,Ui; -yj,Uj) + Bs{-yi,Ui;yj,Uj)) - 

i<j i<j 

^{Bs{xi,ti; -Xi,ti) + Bs{-Xi,ti;xi,ti) + Bs{y„Ui; -yi,Ui) + Bs{-yi,Ui;yi,Ui))] (A9) 

i 

where 

Ba{x,t;y,u) = {C.aL{x,t)C.aL{y,u)) - ^{C,aL{x,t)C,a,L{x,t)) - ^{C.aL{y,u)Ca,L{y,u)). (AlO) 

The computation of this free correlator can also be done in terms of ^-functions, with the result 

B,{x,t;y,u)='-'^^^^^^±^^~^\nFMt-u)+x-y) (All) 



where 



in which a is an implicit cutoff, i.e. z ^ z — ia. The periods of the theta functions are 

'^a = (A13) 

This whole collection of formulas can now be assembled to yield our original correlator. In doing this, phases 
simplify somewhat; in the end, we obtain the general formula 



~.(2n) 

n<jusii i^a. Slim c,a ^"^^f- ■^Za i — r smti ztn 
[Pa+-] ^ [Pa-+] ^ [Pa++Pa-]^ [J l^^a (2a;.)F, (2y,) | ^ (A14) 



Gf^^^y{{x,t;y,u}) = J\f^'^dj2.ai,j:,TiQ{a,T}{{x,t;y,u}) x 

, o n 



a=c,s i=l 

where for ry, = ± we have defined the functions 

P n"j=i Fc{vc{ti - Uj) + rjiXj + ri2yj) 



CV1V2 Yl'l^^.^^ Fc{vc(f.t - tj) + vixt + ■n2Xj)Fc{vc{u^ - Uj) + rjiyi + r)2yj) ' 

-,0"; T- 



n",=i Ps ' ' [vsiU - Uj) + riix, + ri2yj) 



smm Y{'l<j=i Fs'"' {vs{ti - tj) + rjiXi + i]2Xj)Fr''' {vs{ui - Uj) + rjiyi + 772%) 
and where the zero mode part is 



(A15) 



n (U .,\\ - ^2(4 + T,kFR\TMdl\Ts) + Md'c + TckFR\Tc)Mdl\rs) A2dik^R/n ( ..(^^ 

Qi„,r}[ix,t,y,ui) - ^,ir.kpR\Te)d,{0\T.)+MrckFR\re)MO\rs) ^ ' 



in which we have defined = i \ and 



c^c = - ^ XI [e~^^°^'c(ti - Ui) +Xi- yi 



i=l 



(fs = X! '^^"'"sic^iU - TiUi) + GiXi - nyi] . (A17) 



2R . ^ 

1=1 
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One can convince oneself that the correlators obey the correct (anti-)periodicity by shifting imaginary time like in the 
Matsubara formalism, and that moreover they obey the spatial periodicity by shifting distance arguments by twice 
the system size. This will of course be true of all correlators we present, in all differing boundary conditions cases. 

Other correlators, for example of the type where the order of the 'i/^ differs, can be easily obtained from the 
above formulas: the only differences with the formulas that we have presented are some simple signs, in front of 
the coordinates and exponents. We have left these out here to keep the formulas more compact, and since these 
modifications are trivial to work out. 



2. Luttinger liquid, Andreev boundary conditions 



The general correlator is here again defined as 

n n 



(A18) 



i=l 



where we now deal with Andreev boundary conditions, specified by the subindex A. The computation is rather similar 
to the one above; in fact, one can convince oneself rather easily that the general formula becomes a simple variation 
of the one for normal boundary conditions, i.e. 



n 



[Pa-+] ^ 



Gf^ly{{x,t;y,u}) = N'„^6j2.ai,j:,nQ{a,T}{{x,t;y,u}) x 

[^-./^c--]-"^[w.--]^ri§§#MS (A19) 



tJi\Fs{2xi)F,{2yi)\^ 

where the P functions are as for normal boundary conditions, and where the zero mode part is 

n f . - ^2(4 + feXmfc)Mdl\rs) + ^3(4 + fcXmfc)Mdl\rs) .J^w. 

i^i.,r}[ix,t,y,ui) MfcX/2\fc)M0\rs)+MfcX/2\fc)U0\rs) 

in which we have defined fr = i ^"°d and 



(A20) 



dl = -J^Yl [e^^°^c(ti - Ui) +Xi- Vi] . 



(A21) 



The other parameters are defined as in the previous section. 



3. Carbon nanotube, normal boundEiry conditions 

The general correlator for carbon nanotubes in the presence of normal boundary conditions can be written in terms 

of 

n n 

G\Z\i{^'t'yM) ^ iU'^LlM^k) n (A22) 
j=i j=i 

where Sj, fj = 1,2 refer to the nanotube channel index. The computation of this general correlator is very similar to 
the computation presented above for a Luttinger liquid with normal boundary conditions. The final formula is 

g\Z\ ({^' t; y,u})= Af-'SKieinQ\Z\ ii^^ "}) X 

X [P+_]""^ [F+ J""^ [P++^^ct-]"^ [F-+_nVn+-]"'^'ni^c(2x,)Pc(2y,)r'^. (A23) 

i=l 

In the above, dKiein is a shorthand for the Kronecker deltas coming from the Klein factors, meaning that the set of 
(pairs of) indices associated to the ^'^ operators has to be a permutation of the set associated to the ^ operators. 
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For later convenience, wc also define s(si) as +1 for the 1 channel, —1 for the 2 channel (and similarly for s{rj)). We 
have defined the dynamical contributions 



P± - 

ar)ir/2 yi" 



with indices 



icU) = h icU) = itU) = 77 (i) = ^(^jOT-j- 



(A24) 



and = = Va, with a = c, s. The zero mode part is 



(A25) 



(A26) 



where is a four-dimensional ^-function (such multivariable functions are defined in the next appendix). The four- 
dimensional period matrix is (cr^ is the Pauli matrix) 



1 ® (1 + y(1 + e-'«=) - ® (1 + C7')^{1 - e-2«=) + 1 ® (1 - a^)t;. 



and Tc = p . The vector u is 



u = 



f d+ + d- +d+ +dj\ 
dt + d- -dt- d- 
dt-d-+d+- d- 

\d+-d--d++ d- ) 



in which the distances are given by 

n 



Also, we put = = 0, = and fc^ = fc^ = fc^ 



(A27) 



(A28) 



(A29) 



4. Carbon nanotube, Andreev boundary conditions 

For Andreev boundary conditions, we define 



(A30) 



where Sj, rj = 1,2 refer to the nanotube channel index. The computation is very similar to those above, and yields 

[pi_]""^ [p++]"'^ [pi+pi.]""^ [p-+_nvn+-]"'^'ni^c(2a;,)i^c(2yi)i 



sinh 



(A31) 



The zero mode part becomes here 



Q|^;^}({a;,t;2/,u}) 



with period matrix 



n = i 



4R 



1 <^ (1 + + e'«0 - ^1 O (1 + (Ti)y (1 - e2«=) + 1 ® (1 - ai)t;. 



(A32) 



(A33) 
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and Tr = i'^^^ — . The vector u is 



. 1 



/ d+ + d-+d++d- \ 
d+ + d- -d+ -d- 
d+ - d- +d+ -dj 
\d+~d- -dt+dj ) 



(A34) 



The newly defined distances (only d+ is actually needed, as in this particular setup the others coincide with the 
previous ones) are given by 



n 



(A35) 



We have also put xi = X2 = X- 



APPENDIX B: ^-FUNCTIONS, AND LIMITS OF CORRELATION FUNCTIONS 

1. Traditional functions 



Theta functions are defined with the following series: 



^inm -{-i2nu 



iTTTn -\~i2nu 



(Bl) 



The most important property of these functions is their quasi-double- periodicity in the complex plane. We refer the 



reader to the standard literature for a list of all the properties obeyed by these remarkable functions |33 

e and/or zero terr 
periods. We here 

2e-W4 sin(w)-H 



In taking different limits like that of finite size and/or zero temperature, we need to consider various limits of theta 
functions with respect to their arguments and periods. We here collect a few useful formulas: 



9i{u\ia 
9i {u\ia 

92{u\ia 

where ... denote subdominant terms. 



>a^o 2a ^/^e 4a sinh ^—^ 
>„^oo 2e"™/'' cos(u) + ... 



.0 oT^'^e 



1 - 2e"''/" cosh 



1 + 2e-''" cos(2u) 



.0 



1 + 2e"''/" cosh 



2u 



2u 



'a^oc - ' COs(2u) 



1 - 2e- 



.0 a-i/^e-^ 



26-^^/4" cosh (^) 



(B2) 



2. Multivariable theta functions 

For the nanotube, the above ^-functions are not sufficient. It is rather more convenient to use multivariable functions 
531, defined as 



|j-2"j ^ ^ gi7rn*f2n+2m'z 



(B3) 
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where g is an integer (the "dimensionahty" of the function), Q, is a g x g matrix with positive imaginary parts, 
and the sum is taken over all (7-dimensional real vectors with integer coefficients. As for the usual ^-functions, the 
multivariable version obeys two quasiperiodicity properties in the complex plane. These are, for any given m e Z^, 

e{z + Tnn\fi) = e{z\n), 

e{z + TTfimin) = e-''''"'"'"-2''"'''6l(z|n). (B4) 

The limits of low or high temperature are obtained by taking det r2 ^ 00 or det 17 ^ respectively. In the first 
case 

a 

6'(z|f2) l + ^e""-^2cos(2zi) (B5) 

For the high temperature limit we need to use a generalization of the Poisson summation formula to given by 

9{z\il) = , ^ y e-»i(n*-+-*)"-^(M). (B6) 
Vdet(-jn) 

The limit is then easily taken. 

APPENDIX C: NOTES ON THE REAL-TIME PERTURBATIVE FORMALISM IN FINITE-SIZE 

SYSTEMS 

In order to compute thermodynamic physical quantities for a given system at finite temperature, one usually uses 
the Matsubara formalism using imaginary time. If one desires to obtain finitc-tcmpcraturo timo-depcndcnt correlation 
functions, one then has to perform an analytical continuation back to real times, an operation which, done correctly, 
is often inconvenient. 

We wish here to discuss some details of the derivation of the Matsubara formalism from first principles, and show 
that in our geometry of interest, i.e. a one-dimensional finite-size system, there are subtleties one has to bear in mind. 

We adopt the following logic, using a real-time formalism to start with. Let us suppose that at t = —00, the 
Hamiltonian of the system is described by a known operatorial expression Hq. We then adiabatically turn on an 
additional piece Hi with a modulation function f{t). That is, H{t) = Hq + f{t)Hi. The fundamental object that we 
need to consider is the density matrix p{t), which obeys the time-evolution equation 

^^^pit) = [Hit),p{t)]. (CI) 

Defining p{t) = po + pi{t) where po = e~^^°/Zo, we can obtain a self-consistency equation for pi{t): 

pi{t) = -if dti/(ii)e-*^°(*-*i)[iJi,po+Pi(ti)]e*^°(*-*^) (C2) 

J — CO 

which we can solve as a series pi{t) = J2T=i Pi"^(0/^Oj each successive term containing one additional time integral 
and commutator with the unperturbed density operator. 

Let us consider the first correction. For simplicity, we consider a modulation function with periodicity f{t — ifi) = 
f{t), although this is not a fundamental requirement. Using the notation Hi{t) = e^^°*Hie~^^°* , we obtain after 
changing integration variables and using the periodicity of f{t) 



pW(^) = -^e-'3^o 



■ 1—10 M 

i dti- dti 

J — oo—iP J — 



f{h+t)Hi{h). (C3) 



We can then combine the two integrals into one closed contour (using /(— 00) = by construction) and a leftover 
Matsubara piece: 



p^l\t) = 27re-'3^« ^es(/(ti -t- t)Hi{ti) - ie"^^" / dtif{ti + t)Hi{h). (C4) 
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Let us suppose that we are in an infinite (and therefore aperiodic) system. If we choose to turn on our interaction 
Hi at a very remote point in the past as compared to any object we want to calculate (say t = 0), then f{t) will 
have poles only extremely far away, and the first correction will yield expectation values with extremely large time 
difference, which one can safely take to vanish. Moreover, near t — 0, we suppose that f{t) has reached its asymptotic 
behaviour f{t) « 1. The second term then becomes 

P[%tsit) = f drHl{-^T) (C5) 

JO 

which is just the Matsubara contribution. 

This logic, however, fails for a finite-size system described by an initial Hq containing no dissipation. To see this, 
one simply has to realize that in such a finite-size system, correlation functions are by construction periodic in time, 
with a return time of the length of the system divided by the appropriate velocity. For example, one can take as f{t) 
a step function triggering the perturbation Hi at time ^q. In that case, the computation above gives 

p W (t) =e-^"'> (\t [Hi (to - ir) ~ Hi{-it)] (C6) 
Jo 

and no matter how far we put to in the past, the first term never vanishes, since all correlations in the system are 
periodic. In all that we have done, we neglect the contribution from the first term, which we presume "die off" for to 
enough in the past. For the second-order contribution, similar arguments yield 

Plit) = / dTi / dT2f{t - lTi)j(t - IT^)T, [Hi{~lTi)Hi{~lT2)] , (C7) 

^ Jo Jo 

(where Tr is the imaginary-time ordering operator) which we use in the text with the same simplification as used in 
the first-order term, i.e. replace f{t — it) = 1 in the form above. 
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